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ABSTRACT 

The  fractal  dimension  of  a  surface  is  a  measure  of  its  geometric  complexity  and 
can  take  on  any  non-integer  value  between  2  and  3.  Normally,  the  topological 
dimension  of  surfaces  is  2;  however,  their  fractal  dimensions  increase  with  greater 
amounts  of  complexity  or  roughness.  For  example,  a  fractal  dimension  of  2.3  is 
found  to  be  a  common  value  in  describing  the  relief  on  the  earth. 


This  paper  discusses  and  presents  examples  of  an  algorithm  designed  to  measure 
the  fracticality  of  surfaces.  The  algorithm  was  developed  at  The  Ohio  State 
University  and  is  shown  to  be  reliable  and  robust.  It  is  placed  in  an  interactive 
setting  and  is  based  on  the  premise  that  the  complexity  of  Isarithm  lines  may  be 
used  to  approximate  the  complexity  of  a  surface.  The  algorithm  operates  with  the 
following  scenario:  Starting  with  a  matrix  of  Z-heights,  an  isarithm  interval  is 
selected  and  isarithm  lines  are  constructed  on  the  surface.  A  fractal  dimension  is 
computed  for  each  isarithm  line  by  calculating  their  lengths  over  a  number  of 
sampling  intervals.  The  surface's  fractal  dimension  is  the  result  of  averaging  the 
fractal  dimensions  of  all  the  isarithm  lines  and  adding  1. 

Potential  applications  for  this  technique  include  a  new  means  for  data  compression, 
a  quantitative  measure  of  surface  roughness,  and  be  used  for  generalization  and 


INTRODUCTION 


The  problem  of  describing  the  forms  of  curves  and  surfaces  has  vexed  researchers 
over  the  years.  For  example,  a  coastline  is  neither  straight,  nor  circular,  nor 
elliptic  and  therefore  Euclidean  lines  cannot  adequately  describe  most  real  world 
features.  Imagine  attempting  to  describe  very  rough  or  bumpy  terrain  in  terms  of 
classical  geometry.  An  intriguing  concept  proposed  by  Mandelbrot  (1967,  1977, 
1982)  is  to  use  fractals  to  fill  the  void  caused  by  the  absence  of  suitable  geometric 
.  representations.  A  fractal  characterizes  curves  and  surfaces  in  terms  of  their 
complexity  by  treating  dimension  as  a  continuum.  Normally,  dimension  is  an 
integer  number  (l  for  curves,  2  for  areas,  and  3  for  volumes);  however,  fractal 
dimensions  may  vary  anywhere  between  1  and  2  for  a  curve  and  2  and  3  for  a 
surface  depending  upon  the  irregularity  of  the  form.  Although  individual  fractals 
have  been  around  since  the  1900's,  Mandelbrot  was  the  first  to  recognize  their 
applications  outside  of  mathematics. 
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ThU  paper  discusses  an  algorithm,  developed  at  Ohio  State  and  Implemented  In  an 
Interactive  setting,  designed  to  measure  the  fractlcality  of  a  surface.  It  also 
presents  results  from  examining  three  surfaces.' 

DEFINITION  OF  FRACTALS  AND  5ELP-SIMILARITY 

In  Euclidean  geometry  every  curve  has  a  dimension  of  1  and  every  plane  has  a 
dimension  of  2.  This  Is  generally  referred  to  as  the  topological  dimension  (Dt). 
These  dimensions  remain  constant  no  matter  how  complex  or  irregular  a  curve  or 
plane  may  be.  For  example,  the  Rocky  Mountain  area  contains  many  Irregularities, 
but  the  topological  dimension  remains  2. 

In  the  fractal  domain  a  curve's  dimension  may  be  between  1  and  2  according  to  its 
complexity.  The  more  contorted  a  straight  line  becomes,  the  higher  Its  fractal 
dimension.  Similarly,  a  plane's  dimension  may  be  a  non-integer  value  between  2 
and  3.  The  fractal  dimension  for  any  curve  or  surface  is  denoted  by  (D)  and  within 
this  framework!  D  £Dt.  Mandelbtut  (1977)  proposes  the  following  definition  for  a 
fractah  "A  fractal~will  be  defined  as  a  set  for  which  the  Hausdorff-Beslcovitch 
dimension  strictly  exceeds  the  topological  dimension." 

Central  to  the  concept  of  fractals  is  the  notion  of  self-similarity.  Self-similarity 
means  that  for  any  curve  or  surface  a  portion  of  the  curve  or  surface  can  be 
considered  a  reduced  image  of  the  whole.  However,  seldom  In  nature  (crystals  are 
one  exception)  does  self-similarity  occur  and  therefore  a  statistical  form  of  self- 
similarity  is  often  encountered.  In  other  words,  if  a  curve  or  surface  is  examined 
at  any  scale  It  will  resemble  the  whole  in  a  statistical  sense*  therefore,  D  will 
remain  constant.  Brownian  motion  is  an  excellent  example  of  statistical  seif- 
slmilarlty.  Because  of  this  principle,  a  curve  can  be  decomposed  into  N«r 
nonoverlapping  parts  and  each  sutaegment  has  a  length  of  l/r«l/N.  Similarly,  a 
unit  square  can  be  divided  into  N*r  squares,  where  the  similarity  ratio  is  r(N)  »  l/r 
*  1/N  .  In  either  case  the  following  equation  applies* 

I  D-log  N/log  (l/r)  (1) 

and  could  be  called  the  shape's  similarity  dimension.  D  ean  also  be  expressed  as: 

P.log  (N/N0)/Iog(X0/X)  (2) 

where  X  and  X  are  two  sampling  Intervals  and  N  and  N  are  the  number  of  such 
Intervals  contained.  If  a  curve  resembles  a  stralgnt  line  .then  when  the  sampling 
Interval  is  halved,  N  doubles  and  the  proportion  equals  1.  The  majority  of 
cartographic  curves  are  not  straight  lines  and  therefore  N  will  more  than  double 
causing  D  to  be  greater  than  1.  The  principle  of  self-similarity  is  dismissed  by 
Goodchild  (1980),  Hakanson  (1978),  and  Scheldegger  (1970).  Hakanson,  for  example, 
points  out  the  absurdity  of  postulating  the  validity  of  self-similarity  down  to  the 
size  of  the  pebbles  on  the  coastline  and  at  the  molecular  interstices  of  those 
pebbles.  Goodchild  demonstrates  that  although  Richardson  (1961)  found  the  west 
coast  of  Britain  to  have  a  constant  D  of  1.23  over  sampling  intervals  between  10 
and  1000km.,  he  found  the  east  coast  to  vary  between  1.13  and  1.31  for  a  similar 
sampling  interval.  This  suggests  that  whatever  created  the  irregularities  on  the 
coastline  acted  at  specific  scales.  Goodchild  states  that  since  self-similarity  is 
only  one  aspect  of  the  fractal  approach,  it  would  be  unwise  to  reject  the  entire 
concept. 


DEVELOPMENT  OF  THE  FRACTAL  SURFACE  ALGORITHM 

The  following  algorithm  is  based  upon  the  research  performed  by  Goodchild  (1980). 
He  noted  the  earlier  empirical  work  performed  by  Richardson  (I960  and  later 


extended  by  Mandelbrot  (1967).  Rlchardion  measured  the  lengths  o t  several 
frontier*  by  manually  walking  a  pair  of  divider*  along  the  outline  »o  a*  to  count  the 
number  of  step*.  The  opening  of  the  divider*  (n)  was  fixed  in  advance  and  a 
fractional  ildp  wa*  animated  at  the  end  of  the  walk.  The  main  purpose  In  this 
section  of  Richardson's  research  was  to  study  the  broad  variation  of  2n  with  n. 

Richardson  produced  a  scatterplot  in  which  he  plotted  log  total  length  against  log 
step  size  for  five  land  frontiers  and  a  circle.  Mandelbrot  (1967)  discovered  a 
relationship  between  the  slope  (0)  of  the  lines  and  fractal  dimension  (0).  To 
Richardson  the  slope  had  no  theoretical  meaning,  but  to  Mandelbrot  it  could  be 
used  as  an  estimate  of  1-D,  which  leads  toi 


*  D»l-0  (3) 

Goodchlld  computed  the  fractal  dimensions  of  surfaces  by  constructing  contour 
lines  and  calculating  their  lengths.  He  extended  Hakanson's  (1978)  analysis  in  order 
to  calculate  D.  In  computing  length,  Hakanson  counted  the  number  of  intersections 
between  a  grid  and  a  coastline,  and  for  a  complex  curve,  it  is  possible  that  any 
even  number  of  intersections  could  exist.  Goodchlld  points  out  the  scale  of  the 
map,  which  determines  the  size  of  the  grid,  Is  not  strictly  related  to  a  sampling 
Interval.  In  terms  of  D,  it  would  be  more  appropriate  to  count  the  number  of  cells 
.  intersected  instead  of  the  number  of  intersections. 

Goodchlld  estimated  D  for  several  self-similar  surfaces  generated  by  shear 
displacement  (Mandelbrot  1973).  He  selected  a  contour  line  at  a  height  equal  to 
the  mean  of  the  minimum  and  maximum  heights.  Cells  were  aggregated  into  larger 
aggregates  of  3x3,  7x7,  9x9,  and  up  to  19x19,  and  classified  them  as  above  (black) 
orbelow  (white)  the  contour  line.  To  count  the  number  of  boundary  cells  or,  in 
other  words,  compute  the  length  of  the  contour  line,  as  a  function  of  aggregate 
size,  a  count  was  made  of  the  number  of  black  aggregates  containing  at  least  one 
white  cell.  This  is  where  the  cells  are  cut  by  the  boundary.  The  log  of  the  average 
number  of  boundary  cells  was  then  regressed  against  the  lo^of  the  aggregate  size. 
Goodchlld  calculated  0  for  the  contour  line  using  Equation  3  and  extending  that 
same  principle  to  a  surface,  the  following  equation  is  used 

D-2-0  (*) 

The  following  algorithm  developed  by  Shelberg  (1982)  extend*  Goodchild's  work  In 
that  D  can  be  approximated  for  nonseif-similar  surfaces  and  is  placed  in  an 
interactive  setting.  Because  Goodchild's  study  only  dealt  with  self-similar 
surfaces,  he  coula  select  a  contour  line  at  a  height  equal  to  the  mean  of  the 
minimum  and  maximum  heights.  A  self-similar  surface  possesses  the  same 
statistical  properties  at  all  scales  and  at  all  locations;  therefore,  choosing  the 
average  contour  line  is  a  valid  method.  However,  a  surface  may  not  be*  entirely 
self-similar  and  selecting  the  mean  contour  line  may  'produce  a  false  fractal 
dimension.  For  example,  If  a  surface  consisted  of  a  plain  and  a  mountain  range, 
then  selecting  the  average  contour  line  would  neglect  both  points  on  the  plain  and 
mountain  range  and  bias  the  overall  fracticality.  Although  two  of  the  surfaces 
reported  in  this  paper  are  seif-similar,  the  algorithm  Is  designed  to  analyze 
nonself-similar  surfaces.  Instead  of  choosing  only  the  mean  contour  line,  any 
number  of  contour  lines,  up  to  a  maximum  of  200  can  be  used  to  examine  a  surface. 

AUsurfaces  are  in  the  form  of  a  Digital  Elevation  Model  (DEM)  and  the  algorithm 
begins  with  the  user  inputting  whether  the  fractal  dimension  will  be  determined  by 
rows  of  columns.  This  option  is  provided  so  that  a  trend  in  the  surface  may  be 
captured.  The  contour*  interval  is  then  entered  and  the  number  of  contour  lines 
determined.  Next,  the  maximum  cell  size,  which  la  comparable  to  Richardson's 
opening  of  the  pair  of  dividers,  is  input.  The  minimum  allowable  maximum  ceil 


size  U  5  because  any  numbtr  leu  than  5  would  leave  the  linear  regression  open  to 
queitlon. 

For  each  contour  line,  cells  are  aggregated  Into  larger  aggregates,  up  to  the 
maximum  cell  size,  and  clanlfled  ai  above  (black)  or  below  (white)  the  contour 
line.  To  count  the  number  of  boundary  ceils,  ai  a  (unction  of  aggregate  size,  a 
count  It  made  of  the  number  of  black  aggregates  containing  at  least  one  white  ceil. 
Figure  1  represents  a  6x6  surface  In  which  the  contour  tine  of  35  is  denoted  by  the 
darker  line  on  the  surface. 

The  algorithm  starts  with  the  lowest  contour  line,  In  this  case  35,  and  classifies 
each  cell  as  white  (1)  or  black  (2)t  see  Figure  2.  Always  beginning  with  a  cell  size 
of  I,  It  then  compares  each  neighboring  ceil,  along  the  columns,  for  boundary  cells. 
For  example,  the  comparison  between  (1,1)  and  ( l,2)f  and  (1,5)  and  (1,6)  Indicates  a 
boundary  cell  is  not  present.  Conversely,  the  comparison  between  cells  (1,4)  and 
(l,5)j  and  (6,2)  and  (6,3)  indicates  boundary  cells  exist  and  in  these  cases  the 
number  of  boundary  cells  Is  Incremented  by  1.  After  the  entire  surface  is 
examined  using  1  cell  size,  the  search  begins  again  using  a  cell  size  of  2.  Now  cells 
(1,1)  and  (1,3);  and  (6,2)  and  (6,4)  are  compared  and  so  on.  This  continues  until  the 
.  maximum  cell  size  is  reached.  For  each  contour  line  the  same  process  is  repeated. 

(t  is  possible,  for  a  given  cell  size,  that  no  boundary  cells  are  encountered.  In 
Figure  3,  this  occurs  when  the  cell  size  equals  3,  4  or  5.  For  example,  the  only 
possibility  of  encountering  a  boundary  cell  is  in  row  6  when  the  contour  line  equals 
19.  If  the  cell  size  equals  5,  cells  (1,1)  and  (1,6)  are  compared  and  no  boundary 
oalis  are  reported  present.  The  results  remain  the  same  for  the  remaining 
comparisons  within  that  row.  For  any  given  cell  size,  if  no  boundary  cells  are 
encountered,  the  contour  line  is  eliminated  from  the  calculations.  Since  the 
algorithm  regresses  log  (average  number  of  boundary  cells)  against  log  (cell  size), 
without  eliminating  the  contour  line,  but  just  deleting  the  cell  size  from  the 
regression,  could  result  in  1,  2  or  3  points  used  in  defining  the  regression  line  and 
immediately  poses  obvious  difficulties.  Although,  by  eliminating  some  contour 
lines,  the  method  neglects  certain  small  portions  of  the  surface,  it  accounts  for 
greater  variations  than  if  only  the  mean  contour  line  is  selected.  The  strategy  is  to 
choose  the  lowest  contour  interval  with  the  fewest  contour  lines  being  eliminated. 

After  counting  the  number  of  boundary  ceils  per  cell  size  for  each  contour  line,  a 
linear  regression  is  performed.  For  every  contour  line,  log  (average  number  of 
boundary  cells)  Is  regressed  against  log  (ceil  size).  A  0  Is  computed  for  each 
contour  line  and  extended  to  a  surface  by  using  Equation  4.  A  surface's  fractlcality 
Is  the  average  of  all  the  D’s  for  the  Included  contour  lines. 

EXAMPLES  AND  RESULTS 

Two  of  the  three  surfaces  used  as  examples  were  generated  by  the  shear 
displacement  method  and  are  statistically  seii»similar.  Each  is  in  the  form  of  a 
50x50  matrix  and,  on  the  CRT,  are  represented  by  a  perspective  view  with  hidden 
lines  eliminated. 

The  first  surface,  in  Figure  4,  has  a  theoretical  D  of  2.2.  Table  1  shows  a  rather 
large  array  of  D  values  and  the  selection  of  numerous  contour  intervals  is  intended 
to  display  the  stability  of  the  algorithm.  In  an  effort  to  preserve  consistency 
among  the  results  and  not  fall  below  the  minimum  acceptable  maximum  ceil  size, 
in  the  following  examples,  the  maximum  ceil  size  Is  3.  In  all  cases,  D  Is  computed 
along  the  rows. 


Figure  4.  30  x  30  surface  with  a  theoretical  D  of  2.2. 

MlnUsm  and  naxlaua  elevatlona  art  >503  and  269. 

*  i 


Figure  3.  30  x  90  aurfaaa  with  a  theoretical  D  of  2.6. 
Min  1mm  and  naxlaun  elevations  art  *42  and  41. 


MOTKt  Heights  ara  unletsss, 


Contour 

Interval 

Maximum 

Cell  Size 

£ 

No.  of 

included  lines 

No.  of  lines 
not  included 

5 

5 

2.1915 

146 

8 

10 

5 

2.1869 

72 

4 

20 

5 

2.1752 

36 

2 

22 

5 

2.2071 

55 

1 

24 

5 

2.2000 

50 

1 

50 

5 

2.1895 

24 

l 

50 

5 

2.2565 

14 

0 

100 

5 

2.2145 

7 

0 

200 

5 

2.2268 

3 

0 

515 

5 

2.1289 

1 

0 

Tabic  1.  Result!  from  examining  a  50x50  surface  with  a  theoretical  0  of  i. 2 

A  D  of  2.1915  would  be  selected  to  represent  the  surface  because  of  the  large 
number  of  contour  lines  (146)  used  in  the  calculations  versus  only  8  contour  lines 
being  eliminated. 

Figure  5  shows  a  50x50  surface  with  a  theoretical  D  of  2.6.  Table  2  displays  the 
results  over  a  number  of  contour  intervals. 


Contour 

Maximum 

No.  of  - 

No.  of  lines 

Interval 

Cell  Size 

/  e 

2.6720 

included  lines 

not  included 

2 

3 

34 

•  7 

5  - 

5 

2.6671 

14 

2 

7 

3 

2.6679 

10 

l 

10 

3  . 

2.6290 

7 

0 

20 

3 

2.6590 

3 

0 

4L 

3 

2.5985 

'  1 

0 

Table  2.  Results  from  examining  a  50x50  surface  with  a  theoretical  D  of  2.6  ' 

A  D  of  2.6720  would  be  selected  because  54  contour  lines  are  Included  in  the 
calculations  versus  only  7  being  eliminated. 

The  third  surface  is  also  50x50,  and  is  a  Defense  Mapping  Agency's  standard 
elevation  matrix  over  a  small  portion  of  Southern  Nevada.  The  minimum  and 
maximum  heights,  in  meters,  are  2011  and  2576.  Although  the  degree  of  self- 
similarity  Is  not  known,  the  D  values,  in  Table  5,  are  extremely  stable.  A  D  of 
2.1895  would  be  selected  to  best  approximate  the  complexity  of  the  surface.  The 
actual  surface  is  not  shown  because  the  hardcopy  Is  illegible. 


Contour 

Maximum 

No.  of 

No.  of  lines 

Interval 

Cell  Size 

D 

included  lines 

not  included 

5 

5 

2.1893 

66 

6 

15 

5 

2.1905 

22 

l 

20 

5 

2.1707 

16 

1 

150 

5 

2.1701 

1 

0 

Table  5.  Results  from  examining  a  "real-world"  50x50  surface. 
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flgurt  1.  6x6  aurface  with  a  contour  Una  of  35. 
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figure  2.  6x6  aurfaca  vhara  balov  contour  Una  “  X 

and  abov#  contour  llna  *  2, 
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Flaura  3.  6x6  aurfaea  with  a  contour  Una  of  19. 
A  eaaa  whara  no  boundary  ealla  art  anco^o'  ad  for 
call  alcaa  3,  4  and  5» 
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*.  SUMMARY  AND  CONCLUSIONS 

The  results  demonstrate  the  stability  of  the  algorithm  In  which,  over  a  number  of 
contour  or  sampling  intervals,  the  fractal  dimension  remains  reasonably  constant. 
Other  30x30  self-similar  -.urfaccs  which  were  examined,  ranging  In  D  from  2.1  to 
2.9,  produced  similar  results.  By  using  a  number  of  contour  lines,  the  maximum 
variation  In  the  complexity  of  the  surface  can  be  captured  and  therefore  D  can  be 
closely  approximated. 

Potential  applications  of  this  algorithm  Include  a  new  means  to  compute  surface 
roughness,  a  means  to  measure  the  amount  of  surface  generalization  or  filtering 
attained,  and  a  method  to  store  a  compressed  surface  so  that  greater  complexity 
can  be  Introduced  to  match  a  target  surface  with  a  specific  fractal  dimension. 
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